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Quantum fluctuations in an anharmonic superconducting circuit allow for frequency conversion of 
individual incoming photons. This effect, linear in the photon beam intensity, leads to ramifications 
for the standard input-output circuit theory. We consider an extreme case of anharmonicity in 
which photons scatter off a small set of weak links within a Josephson junction array. We show 
that this quantum impurity displays Kondo physics, and evaluate the inelastic photon scattering 
cross-sections. These cross-sections reveal many-body properties of the Kondo problem, which are 
hard to access in its traditional fermionic version. 

PACS numbers: 74.81. Fa, 72.10.Fk 



Propagation of small-amplitude electromagnetic waves 
through an optical system or a passive microwave circuit 
is conventionally described in terms of transmission and 
reflection amplitudes, or, equivalently, complex admit- 
tances. Considered classically, the wave propagation can 
be calculated using input-output theory [3, 0]. In the 
absence of dissipation, the transmission t{uo) and reflec- 
tion r{u) amplitudes for a photon of frequency uj satisfy 
the unitarity condition, |i(w)| 2 + |r(cj)| 2 = 1. It is of- 
ten tacitly assumed that this description applies in the 
quantum limit too. While this is indeed true if the cir- 
cuit is harmonic, the presence of anharmonic elements 
modifies the picture qualitatively: a photon of energy fko 
may "split" into several ones of smaller energy; unitarity 
is violated in the elastic channel, |i(w)| 2 + |r(w)| 2 < 1. 
The photon frequency conversion results in a finite dissi- 
pative part of the admittances despite the system being 
free of dissipative elements. These features appear in a 
quantum circuit containing even a single or a small group 
of anharmonic elements, a "quantum impurity" . 

In this paper, we consider propagation of microwave 
photons (oscillations of charge and superconducting 
phase) along an array of Josephson junctions interrupted 
by a capacitive element, see Fig.Q] If Josephson energies 
were all large with respect to charging energies for each 
of the tunnel junctions, the system would be effectively 
harmonic, and photon scattering off the central capaci- 
tive link would be purely elastic. We will rather assume 
the Josephson energy to be large for all the junctions ex- 
cept for the two closest to the capacitive link. These two 
junctions, together with the two superconducting islands 
they single out, form a quantum impurity which causes 
inelastic photon scattering. The quantum impurity is of 
the Kondo variety [3[ , where the two values of the polar- 
ization charge of the said two islands play the role of the 
Kondo spin. However, photon scattering is quite different 
from the electron scattering in the conventional Kondo 
problem Q. We find that the photon elastic transmis- 
sion and reflection coefficients, as well as the total inelas- 
tic scattering cross-section ^(uj), are related to the local 
"spin" susceptibility XzzW). We then study the spec- 
trum 7(u/|cli) of photons at frequency uj' generated by 




FIG. 1. Diagram of the system. The dotted box denotes the 
quantum impurity. See the text for further details. 



inelastic processes from incoming photons at frequency 
uj. The spectrum peaks as a function of uj' at the Kondo 
energy scale. At uj — uj' <C T k or uj' <C T K the behavior 
of j(uj'\uj) provides a direct access to corrections to the 
Nozieres fixed-point Hamiltonian. We provide technical 
details in the Supplemental Material (SM) [f|. 

Assuming that the superconducting gap is larger than 
any other energy scale, the only relevant degrees of free- 
dom are the number of Cooper pairs rii on island i and 
the corresponding superconducting phase (fi, obeying 
[ipi, Uj] = iSij. The array Hamiltonian is 

H = J2 [2e 2 (m - n°) (CT 1 ) . . (n, - n°) - E*j cosfo - Vi ) 
i,j 

(1) 

where E j and are the matrices of Josephson couplings 
and capacitances, respectively. We will assume nearest- 
neighbor Josephson couplings, and ground- and nearest- 
neighbor capacitances, whose values can be inferred from 
Fig. Q] The gate-induced charge offset on the z-th island 
is n® = Cf Vf/(2e) with Vf and Cj being the gate voltage 
and capacitance to the ground, respectively. 

Away from the quantum impurity the array is uniform: 
except for the quantum impurity islands, all Josephson 
couplings are Ej, and all capacitances to the ground and 
junction capacitances are C g and C, respectively. Prop- 
erties of the uniform array are controlled by two ratios, 
Ej/Ec g and Ej/Ec, of Ej and two charging energies, 
E c = (2e) 2 /(2C) and E c = (2e) 2 /(2C 9 ). Typically 
C/C g > 1 (it is - 10 2 in That allows having the 



2 



impedance of the array Z = [H/ (2e) 2 ] yj2Ec g /Ej on the 
order of the resistance quantum Rq — irh/(2e 2 ), while 
keeping the amplitude of phase slips A ~ e ^V 32E j/ E c 
exponentially small [6(. In an array of length L < a/ A 
(a is the array spacing) the Josephson energy can be re- 
placed by a quadratic term. In addition, in the long 
wavelength limit we may use a continuum description for 
the array [?J (except for the impurity) in terms of Bose 
fields 4>i(x) and pi(x) which represent, respectively, the 
superconducting phase (whose gradient is proportional 
to the electric current) and charge density (in units of 
— 2e per period of the array) in lead I — L, R, obeying 
[4>e(x), p£>(x')] = iSwS(x - x'), 



leads - ^2 „ 

e=L,R 







<[d x 4>i{x)] 2 + - [irp e (x)} 2 \ dx. (2) 
.9 J 



The array is characterized by the velocity of plasmons 
v = a,y/2EjEc g , and by g — Rq/(2Z). C does not af- 
fect excitations of wavelength well exceeding a^/C jC g . 
Thus, the linear dispersion waveguide Hamiltonian ([2]) 
is limited to frequencies within a bandwidth ujq ~ 
(v/a)y/CJC (See SM, Sec. SM.A @). 

Let us now turn to the quantum impurity, islands 
L and R in the dotted box in Fig. [T] We derive its 
low-energy Hamiltonian under the realistic assumptions 

c ■ c a L ,c R 
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C g and C L ,C R ~ yJCCg (See 
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SM, Sec. SM.A @). When the charging energy £™ p = 
(2e) 2 /[2(C L + C R )] [here l/C e = l/C e + is 
large with respect to the Josephson energies Ej' n , the 
total impurity charge «l + n R is quantized. If the gate 
voltages are set to {C 9 L V[ + C R V R r )/(2e) = 1, then to 

L R 

lowest order in Ej ' the islands are restricted to the two 
charging states \0lAr) ( n L — 0, n R = 1) and |1l,0.r) 
(til = 1, n R = 0). We label these two configurations by 
the states of a pseudo-spin, S z = {til — n R )/2 = ±1/2, 
so that S+ = \1 l ,0r){0l,1r\, S- = (S+)\ obeying 

L R 

the standard commutation relations. Finite E 'y allows 
for switching between these two states through virtual 
states with energies of order _E™ P . Eliminating these 
by a Schrieffer- Wolff transformation leads to an effective 
low-energy Hamiltonian (See SM, Sec. SM.A @) 
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and \ LR = C L C R /[(C L + C R )C LR ] ~ ^/UJC « 1. The 
first term in Eq. ([3]) accounts for flips of the pseudo-spin, 
which are accompanied by transfers of discrete charge 



±2e between the two leads [8]. The second term is a 
capacitive coupling between the impurity and the leads. 
The third represents the effect of a gate voltage bias be- 
tween the impurity islands. Hamiltonian ([3]) clearly in- 
troduces anharmonicity into the system. 

Applying the transformation H — > U^HU with U = 
e -%l<l>L(o)-4>R(o)]S z ^ ^.jjg Hamiltonian acquires the form of 
the spin-boson model with Ohmic dissipation [t| [Toj : 



X—c.s 



dx4>\{x) + [irpx(x)} 2 > dx 



B z S z -E L j K S x -Ttvap s {0)S z , (5) 



where p s (x) = [a L p L {x) - a R p R {x)]/(a^/g) and tj> s (x) = 
\f§[ a L4>L{x) — a R 4> R (x)]/a are, respectively, the "spin 
density" and its canonically conjugate momentum field. 
The "charge density" and its conjugate field, p c {x) = 
la R pL(x) + a L p R (x)]/(ay/g) and (f> c (x) = y/g[a R 4> L (x) + 
ttL^j;(a;)]/a, decouple from the impurity spin. The pa- 
rameters aj,. R a nd the coupling parameter a in Eq. ([5]) 
are given by [111 ] 
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(6) 



The spin-boson Hamiltonian ([Q is equivalent 0, [l(| to 
the single-channel Kondo model [3J , describing a localized 
spin exchange-coupled to a bath of noninteracting spin- 
1/2 fermions with bandwidth uiq, 



k % a— "f ,4- k,<r,k' ,ct' 

+ W S - E cL<*>Ck'^+K.c-B z S z , (7) 



where t* a , are the Pauli matrices, I z — 27ru(l — a/y/2), 
and I xy = 2-KaEj R . Given the smallness of Ej R [cf. 
Eq. (0}], isotropic exchange (I xy = I z ) corresponds to 
a 2 m 2 (i.e., g w 1, since Xl R <C 1). The standard 
Toulouse point, where the Kondoproblem is equivalent to 
a noninteracting resonant level occurs at a = 1 

(i.e., g w 2); this point of strongly anisotropic exchange is 
hardly accessible in an electronic realization of the Kondo 
model. In all these cases the Kondo couplings still flow to 
the same strong-coupling fixed point as in the standard 
isotropic case. 

The Kondo impurity is locked into a singlet with its 
environment at energies below the Kondo temperature 
Tk- We define it through the inverse static local impurity 
susceptibility, T^ 1 = d{S z )/dB z \ 

b z =t=o- To the leading 
order in I xy oc EyE R it is given by [l2[ 



T K = c(a)w 



-1 



2/[2-a 2 
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c(a) ~ 1, (8) 



with c(0) = 1. For the strong-coupling physics to show 
up the leads should be longer than v/Tk |13j . 
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We now examine the ac transport properties of the cir- 
cuit. The quantum impurity causes elastic and inelastic 
scattering of incoming microwave photons. The former 
is characterized by the elastic T-matrix Tjh^ui), defined 
as usual by the relation between the single photon prop- 
agators in the presence and absence of the impurity (see 
SM, Sec. SM.B @). It has the structure 



- Uuj) = 



ri(w) - 
ti(w) 



tr{oj) - 



1 



(9) 



with tg(ui) [rt(u>)] being the transmission [reflection] am- 
plitudes for a photon of frequency uj incoming in lead 
I. There is another way to look at elastic scattering, 
which is equivalent to the previous one at zero tempera- 
ture and can serve as its generalization to nonzero tem- 
peratures: one may add to the Hamiltonian ©-© the 
term H ac = 2Vq cos(oji)pi(xi n ), describing an ac gate 
voltage coupled to the island at in the left lead. This 
perturbation generates waves propagating to the left and 
to the right. The latter will scatter off the quantum im- 
purity. The transmission amplitude ti,{bS) [14| for waves 
coming from the left is then the ratio of the transmitted 
and incoming average currents, 



x 4>R (^out ) j Pl{ x 'u 



(or 



(10) 



where transmitted current is measured at x ou t in the 
right lead. In the second equality we have written the 
transmission coefficient as the ratio of the conductance 
of the system with impurity (double angular brackets de- 
note retarded correlation functions) and the correspond- 
ing quantity for a uniform array. Thus, |ti(w)| < 1 and 
is independent of Xin i0 ut , as required. We can write down 
similar expressions for the other scattering amplitudes. 

The equations of motion for the bosonic operators al- 
low us to derive a relation 



T^H - (-l) s <>«- 1 wauat>Xz*(u), 



(11) 



between all the elements of the elastic f matrix and 
the local dynamic differential spin susceptibility of the 
Kondo problem ([7]), Xzz{u) — ((S z ; Szfj^. Thus, a simple 
ac transport measurement on this system yields the dy- 
namic susceptibility of the Kondo model, which is hard 
to access in the electronic realizations of the Kondo ef- 
fect: in those systems charge transport is weakly-coupled 
to the spin dynamics, whereas in our system S z is actu- 
ally the electric polarization of the quantum impurity. 
An incoming electromagnetic wave will generate an ac 
voltage difference ( "magnetic field" ) on the "spin" . The 
impurity electric polarization will oscillate in response 
[through Xzz{u)} and emit the scattered waves. 

The frequency dependence of \zz at low temperatures 
(T <C Tk) is non- monotonia We will concentrate on 
small "magnetic fields" [cf. Eq. flU], B z -C T K , where 
Kondo physics is most clearly manifested. The imag- 
inary part of Xzz{u) has a maximum at w ~ Tk while 
Re[xzz(w)] alternates its sign. These features sharpen up 



at smaller a into a narrow resonance of width ~ o?Tk 
at a -C 1 [10]. At low frequency u) <C Tk and arbitrary 
a the susceptibility approaches a real constant, 



Xzz(uj) = xo ( «, 



El 

T k 



+ nra u>xo ce 



El 

T k 



, (12) 



where xo(<*, B z /Tk) = d(S z )/dB z is the static local dif- 
ferential susceptibility, with Xo( a i 0) = l/T K - The coeffi- 
cient of the dissipative, linear- in- frequency term, is fixed 
by the Shiba relation [SEl (See SM, Sec. SM.C @). At 
high frequencies, uj 3> Tk,B z , we can use perturbation 
theory in I xy cx E^Ef to find [IB] 



Xzz{uj) = i 



.7r f{a) ( T K 



2-a A 



a>l, (13) 



where /(a) = -2sin(7ra 2 /2)r(l-a 2 )/{7r[c(a)] 2 - Q2 }. At 
a < 1 the imaginary part of Eq. (fTB"]) still describes 
Im[xz2(w)], while the real part is dominated by an- 
other term, Ke[x zz (uj)} = — g(a, B z /Tk )Tk/uj 2 , with 
g(a,B z /T K ) ~ 1 for B z < T K . At w » T K ,B Z , the 
photon reflection coefficient |r^(w)| 2 in the elastic channel 
approaches 1, while the transmission coefficient |^(oj)| 2 
scales as (Tk/uj) 2 ^- 2-01 ' for a > 1 and as q, a {Tk/uj) 2 for 
a < 1. 

Let us now turn to inelastic photon scattering. Using 
Eq. (fTTj) , the total probability of an incoming photon to 
be scattered inelastically is 



7« (a>) = 1 - \ri(w)\ 2 - \ti(uj)\ 2 

= 47ra 2 cjlm [x zz {uj)\ - 47r 2 a 2 a 2 w 2 \x zz (uj)\ 



(14) 



This quantity would be zero for a harmonic system, but 
is nonzero in general (See SM, Sec. SM.C [H]). As a 
matter of fact, for uj 3> Tk we may use Eq. (| 13[) to find 
7^(cj) ~ a 4 (Tx/w) 2 ~ Q , which is parametrically larger 
than the elastic transmission coefficient |i^(a;)| 2 for any 
a. 

The measurable characteristic of the inelastic processes 
is the spectrum of emitted photons 7^/u(o/|o;), where 
7^/|^(a/|w)da/ is the average number of photons in the 
frequency interval [a/, a/ + da/] emitted into lead £' per 
each incoming photon at frequency uj in lead I (assuming 
the incoming intensity is weak enough that processes in- 
volving two or more incoming photons can be neglected) . 
This quantity is a sum over the cross sections of all the 
possible multiphoton inelastic processes where one of the 
outgoing photons has frequency a/, while integrating over 
all the other outgoing photons. It can also be related to 
local impurity correlators (See SM, Sec. SM.D [5]). En- 
ergy conservation leads to the relation 



oo 

^ J uj'j e <\e(uj'\uj)duj' = ojje(uj). 



(15) 



e'=L,R 



For uj,uj',uj — uj' 3> B z .Tk the spectrum j£tu(uj'\uj) 
can be found perturbatively in I xy cx EjEf (See SM, 
Sec. SM.E 0]), 
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lt'\t{u'\u) 



Luuj' \4-7ra / 



Im [x+_(w — uj')] (o(uj - w'){[l + ^s(w')] [1 + nfl(w — a/)] — n B (uj')n B (ui - u/) j 
+0(w' - w)|n B (w') [1 + - w)] - [1 + n B {uj' - cj) jj 

+Im[x+_(u; + a/)] {«b(w + w') [1 + ns(w')] - [1 + n fl (cJ + u/)] - ( 16 ) 



where ns(u) = l/(e w / T — 1) is the Bose distribution, and 
x'lLH = | e W s (0) ;e -i^(0)| ) calculated for vanish- 
ing coupling to the impurity. The different terms in this 
equation account for all the possible multiphoton scat- 
tering processes. For example, the first term on the first 
line describes a process where an incoming photon at 
frequency ui is absorbed by the quantum impurity, which 
in turn emits a photon at frequency uj' < uj [hence the 
stimulated emission factor 1 + n B {uj')], plus additional 
photons whose energies sum up to uj — w'. It can be 
shown that the factors depending onu-w' can be writ- 
ten as the sum over the probabilities of distributing the 
energy uj — ui' among any number of photons (See SM, 
Sec. SM.E [|). At T = Eq. CE6|) yields (for uj' < u) 



lt<\t{u'\u) = 7T a t a t ,f(a) 



T K 



, (17) 



with f(a) = sin[7r(a 2 — l)/2]/(a). This result, together 
with Eqs. (|T5 |) -([T1 ]) . obeys the sum rule §T5§ to the lead- 
ing order in Tk 

If any of the energies uj, uj' , or uj — ui' becomes less 
than Tk, perturbation theory in I xy is no longer valid. 
To derive the behavior of 7^/1^(0/ |w) in these regimes, 
let us start from the case when all the frequencies are 
small, and the dynamics is governed by the strong cou- 
pling fixed point. At low energies the impurity is screened 
and disappears from the problem. According to the 
Nozieres Fermi-liquid description, it leaves behind (at 
B z — 0) local scattering potential and interaction be- 
tween the fermions of Eq. ([7]), mediated by virtual fluc- 
tuations of the Kondo impurity. Upon bosonization, 
the leads are described by the first term of Eq. ([5]) 
while the local potential and interaction acquire the form 
H sc = h(a)v 2 p 2 s (0) /T K This is the lowest order 

term allowed by symmetries; for example, the spin den- 
sity oc /5 S (0) cannot appear in odd powers due to the 
time reversal symmetry of the Kondo model, represent- 
ing the equivalence of the two impurity states in Eq. ([5]) 
at B z — 0. The term H sc is harmonic; in order to 
study inelastic effects one needs to consider higher-order 
terms. In the absence of magnetic field, a quartic, four- 
photon term H4 = Z 4 (a)w 4 pf (0)/T K is the lowest an- 
harmonic term allowed, while with magnetic field three 
boson scattering H 3 = l 3 (a)B z v 3 p^{0)/T^ is possible 
[lp(a) ~ a p for a <C !]• Fermi's golden rule then leads to 



(for J < lu < T K ) 



= of of, 



a B {a)Bl 



j (a) (u- 



T* 



(18) 

(the coefficients 0,3, u(a) are evaluated in the SM, 
Sec. SM.F Q for small a). The linear in u>' behavior 
of the spectrum for small uj' stems from the expansion 
coefficients of p s (0) in terms of the boson operators. The 
same applies for small lu — uj' if B z 7^ 0, since then all that 
energy may be carried by a single photon; when B z = 
the energy uj—uj' has to be distributed among at least two 
bosons, leading to a higher order behavior ~ (uj — uj') 3 . 

Returning to the high frequency regime uj ^> Tk, the 
behavior near the edges of the spectrum in uj' is the same 
as for uj <C Tk, since at uj ~ Tk a crossover, rather 
than a singularity, occurs. Thus, while Eq. (|T7| applies 
as long as both ut'.ui — uj' ^> Tk, for small uj' one has 
j£iu(uj'\u)) oc uj', whereas for small uj — uj' > 



■j e , le (uj'\uj) 



2 2 
OLdOLd, 



(uj-u,') a' B (a)B 2 z +aUa)(u-u'y 



uj*T* k 



(19) 

(See SM, Sec. SM.F @, for the small a values ofa' B bJ (a)). 
The resulting behavior is depicted in Fig. [2] 

To conclude, we have considered the scattering of mi- 
crowave photons propagating along an array of supercon- 
ducting islands by a localized anharmonicity. We have 
shown that, contrary to the assumptions of input-output 
theory, linear response is typically dissipative, and inelas- 
tic scattering is therefore significant. The photon scat- 
tering provides direct access to the dynamics of quantum 




FIG. 2. Qualitative behavior of the inelastic spectrum for 
u> 2> Tk- The peak at u>' ~ Tk becomes sharper, with width 
~ a 2 T K , for a < 1 (See SM, Sec. SM.F @|). The peak at 
uj — uj' ~ Tk disappears for a > 1 (cf. Eq. (|17[) ; compare SM, 
Figs. S2-S3 H). See the text for further details. 
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impurity. While we have concentrated on a Kondo sys- 
tem, these conclusions should apply to other types of 
quantum impurities. Finally we note that elastic scatter- 
ing in this system in the limit a«l was recently studied 
in Rcf. 18. Inelastic scattering, whose probability is small 
in that limit (See SM, Sec. SM.F [5|), was ignored there. 
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SUPPLEMENTAL MATERIAL 



In the Supplemental Material we go into some technical details of the calculations, which were omitted in the main 
text. 



SM.A. DERIVATION OF THE EFFECTIVE LOW-ENERGY HAMILTONIAN 



In this Section we will outline the derivation and range of validity of the effective impurity Hamiltonian, Eqs. ([2])- 
(j4]) in the main text, starting from the general array Hamiltonian, Eq. ([!]) in the main text, with the parameters 
given in Fig. [T] in the main text. While this can be done in the general case, the resulting expressions would be 
quite cumbersome. Therefore, we will concentrate on the typical regime of parameters for realistic systems [Sl| . In 
particular, inter-island capacitances are typically much larger than the capacitances to the ground, and the impurity- 
lead capacitances are smaller than other inter-island capacitances: C ~ Clr ^ Cl ~ Cr ^S> C g ~ C£ ~ C g R . As we 
will see in the following [cf. the discussion after Eq. (|S12|) ]. the optimal value of the impurity-lead capacitances Cl,r 
is of order ^/CC g , which we shall assume. In the following we will only keep terms to the lowest nonvanishing order 
in the corresponding small ratios. 



1. The inverse capacitance matrix 



In order to write down the explicit form of the Hamiltonian, Eq. (fTJ) in the main text, one needs to invert the 
capacitance matrix C mm >, where the capacitances can be read off from Fig. Q] in the main text. This can be done 
similarly to the calculation of the Green functions of a noninteracting tight binding model [S2j , where the capacitances 
to the ground take the place of the onsite energies, and the inter-island capacitances are analogous to the hopping 
matrix elements. The presence of large inter-island capacitances makes the inverse capacitance matrix long-ranged: 
for a uniform lead one has [k = 1 + C g j (2C)] 



n(o) 

mm' 



2CV> 



(k — \/ r 2 — l) 

K 2 — 1 V ' 



\m— m'\ 



2^CC g 




1 - \ —^r 



(SI) 



whereas for a half-infinite lead, 



hi 




(k - Vr 2 - l) ' m m ' + (k - \/r 2 - 1) 



I m-\-m' — 1 1 





m+m — 1| 



(S2) 



with m, m' > 0. 

We can now write down the elements of the inverse capacitance matrix of the system in the presence of the quantum 
impurity, which appear in the first term of the Hamiltonian, Eq. ([T]) in the main text. The impurity sub-block of the 
inverse capacitance matrix is given by 













LL 




LR ] 


c- 1 ; 




c- 1 " 




RL 




rrJ 



'ci + c L + c LR -{c L f [c- 1 ]^ 

v —Clr 

i A i\ i 



C h 



C L + C R 



Clr{C l + C R ) 2 



-C LR ^ 
Clr -(Cr) 2 [C- 1 ]^ 

C l Cr\ 
ClCr C 2 )' 



C R 



where Ce = Ci^JCC g j (Ci + y/CC g ), and where 

r c -n*M = r c _ii 

L J mm' L J 



hi 

mm' 



hi 

71 1 1 



C t 



■ 

L J lm' 



1 + Ci [c-^\ 



(S3) 



(S4) 



are the elements of the inverse capacitance matrix of a half-infinite lead terminated by a capacitance to the ground 
whose magnitude is C(. The first term in the last line of Eq. (|S3|) dominates the dynamics of the impurity total 
charge, while the second governs the behavior of its polarization. 
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The impurity-leads elements of the inverse capacitance matrix appearing in Eq. ([T]) in the main text are given by 

whereas the lead-lead elements are modified to (in the following subsection we treat the lead dynamics in a Lagrangian 
formulation, and thus do not use this formula; it is given here for reference): 

[C- 1 ]^, = [C-^i, + [C-j^QCC-']-;,^ [C-]^ . (S6) 



2. Validity of the low-energy Hamiltonian of the leads 

As Eqs. (|Sl|) -(|S2 j) show, the large inter-island capacitances result in a long range of the inverse capacitance matrix. 
For a uniform array this may be ignored as long as one is interested in modes with wavelengths longer than the charge 
screening length a^/C/Cg [see also Eq. (jSlOp below], leading to the low-energy effective leads Hamiltonian, Eq. ([2]) 
in the main text. However, the situation is more complicated in the presence of the nonuniformity created by the 
quantum impurity. The capacitive coupling to the impurity modifies the dynamics of the leads electromagnetic modes, 
allowing for their scattering and transmission between left and right even for Ej = Ef = 0, i.e., in the absence of 
the quantum impurity dynamics. In this subsection we will show that these effects can still be ignored, and Eq. ^fy 
in the main text may still be used, at energies lower than wo ~ {v/a)yjC g /C. 

Let us therefore examine the case Ej = E R — 0, assuming further right- left symmetry C£ = C 9 R — Cq, Cl = Cr = 
Co (effects of right-left asymmetry will be discussed below). After replacing the Josephson couplings by quadratic 
terms, as appropriate for EjC/(2e) 2 3> 1, the symmetric and antisymmetric modes with respect to the center of 
the array decouple. Relabeling the islands to the right/left of the impurity by m = ±1,±2, respectively, the 
symmetric and antisymmetric modes are defined by 



V2 



m > 



v/2 



(S7) 
(S8) 



and similarly for the operators n^. Their dynamics is governed by the Lagrangian 



m>0 





(S9) 



where Cq + = Cq, Cq = 2Clr + Cq. The eigenfrequencies of the system are then 

l(2e) 2 Ej sin(fca/2) 



while the eigenmode expansion is: 



u{k) 



1 + 4^sin 2 (fca/2) 



(S10) 



k 



ka\ m 

2 



m > 0, 



(Sll) 



where L is the leads length (the allowed values of k depend on the exact boundary conditions at the far end of the lead, 
but this is immaterial for the quantum impurity dynamics we are after), and the scattering phase of the eigenmodes 
is given by 



a 2i{Sf-ka) _ 



g± -ika/2 



2 [(2e) 2 Ej - uj 2 (k)C] sin(fca/2) - iw 2 {k)C 9 1 ± e 
2 [(2e) 2 Ej - w 2 (fc)C] sin(fca/2) + iw 2 (fc)Cf ± e ifca / 2 



where Cf ± = C g + Cq ± Co/(Cq ± + Co) is the effective total ground capacitance of the island 
\ S k\ > \ 5 k\- The P hase shifts are negligible for ka < (2e) 2 Ej/(v 2 CD - C g /C - 



(S12) 
1. Therefore, 



8 



Therefore, the use of the low energy effective Hamiltonian, Eqs. (HI)-© in the main text, as well as ne- 
glecting of scattering of photons by the impurity capacitances, are justified only at frequencies smaller than 
m&x[(v/a)y/C g /C, (v/a)C g /Co\. Choosing Co of the order of yJCC g is optimal, as mentioned above, in the sense of 
matching the two cutoffs and thus not "wasting" frequency range. In this low frequency limit the eigenmode expansion 
of the occupancies ra^j — DCq / drf>^ takes the form 



± _ /2 



A,- 



ka 



( I s 

m 

V 2, 



m > 1, 



C„ 



2 

J^2 n k cos(fea/2), 

k 



(S13) 
(S14) 



i.e., only the behavior at m = 1 is significantly affected by the inter- island capacitances. If we lift 
the restriction of right-left symmetry, a similar calculation shows that at frequencies much smaller than 
max[(v / a) y/Cg /C , (v/a)C g /CL r] all the above essentially remains the same, except that Cq is replaced by 
2C L C R /(C L + C R ) [cf. Eq. JS15]) below]. 



3. The quantum impurity Hamiltonian 



Building on the basis laid down in the previous subsections, we will now write down the quantum impurity part 
of the Hamiltonian, Eq. (fTJ) in the main text, at frequencies smaller than ojq ~ (v/a)y/C g /C. Let us start from the 
charging part. The inter-impurity capacitive coupling is given by Eq. (lS3|) . Using Eqs. (jS2j) - (jS5l ) and (jS13|) (|S14[) . 
the effective impurity-lead capacitive coupling at low frequencies assumes the form 



(2ef 
V2 



£ 

e,m>o 

2 



civ? 



2e 



| 2C L C R \ 
c g c LR {c L + c R y V" (C L + C R ) y/CC~ g ) 



(2e) 2 C L C R 
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C L [n L - 



cm 

2e 



Cr n R 



LmR) n ™ I 



C 9 V 9 
R R 



2e 



where 



Plm( x ) = -j=^ J2( n k T n k ) cos(A;a;), 



[p L (0) - P r(0)} (S15) 



(S16) 



are the fields occurring in the continuum version of the lead Hamiltonian, Eq. © in the main text. 

Combining Eqs. (|S3I) and (|S 15[) together with the impurity-lead Josephson coupling, the impurity Hamiltonian 
assumes the form 



-"imp — 



(2e) 2 



n L + n R 



c 9 L vl + c R v£ 



2c 



(2e) 2 



2ClC r 



2{C l + Cr) 

| (2e) 2 C L C R 

c g c LR (c L + c R y \ L T (c L + c R )y/cc~ gj 

- E h j cos [<p L - 0i(O)] - Sfcos [ifR - 4> R {0)] . 



2C L r(C l + CrY 



C R n L - C L n R - 



C R CiV[ - C L C* R Vi 



C L [n L - 



2e 



C R I n R 



2e 



2c 

[ P l(0) - PR (0)} 
(S17) 



Here til, '-Pl and ur, (p R are the number and phase operators of the islands L and R, respectively. 



^. Derivation and validity of the effective spin impurity Hamiltonian HJ\) 

We will now outline how the effective spin impurity Hamiltonian, Eqs. ©-(jl]) in the main text, can be derived 
from the more general form (1S17|) under suitable conditions. 

As mentioned in the main text, the quantum dynamics of phases <Pl,r strongly depends on the ratio of the Josephson 
energies Ej' R to the charging energy i?™ p = (2e) 2 /[2(Cl + Cr)]. If the latter is small, phase fluctuations are small 
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and one may expand the Josephson energy part of -ffj mp to second order in the respective arguments. The resulting 
harmonic version of -ffj mp would lead to elastic photon scattering only. In the opposite limit, E 1 ™ 9 ^> Ey , the total 
charge — 2e(n,L + ur) of the two islands is fixed by the large Coulomb energy penalty. If the gate voltages are tuned to 
a total charge of a single Cooper pair, (C£V^ + C R V R ) / (2e) = 1, then til + n R = 1. When Ey are zero, the charge 
of each of the islands can only take the integer values or 1, and does not vary in time. The possible occupancy states 
are thus |0l, 1r) (i.e., Ul = 0, n R = 1) and \1l,0r) (i.e., til = 1, ur — 0). We label these two charge configurations 
by the states of a pseudo-spin, S z = (ul — tir)/2 = ±1/2, so that 5+ = 0_r)(0l, 1r\, S- = (S + Y . Hence, S z and 
S± obey the standard spin commutation relations. 

Finite Ej and Ef allow for switching between the configurations |0l, 1_r) and |1l,0^) (i.e., flipping of the pseu- 
dospin) by virtual transitions to states with til + n R 7^ 1, with energies higher by ~ i?™ p (which is of the order of 
ujq ~ (v/a)^jC g /C for Ej > Ec )■ The charge part of the Hamiltonian (IS17[) can be projected into the low energy 
sector by substituting til r = l/2±S z , yielding the last two terms of Eq. ([3]) . Terms that do not involve the impurity 
degrees of freedom can be gauged out up to a renormalization of the magnetic field B z [corresponding to the term 
proportional to \ 2 LR in Eq. (@} of the main text] . 

As for the Josephson part of Eq. (|S17[) . one may perform a Schrieffer- Wolff transformation [S3j in order to account 
for processes involving high-energy virtual states. This results in the first term of Eq. ([3]). Here it should be 
noted that the Schrieffer- Wolff transformation also yields terms containing S z pl,r(0)- These would have amplitudes 
~ (Ej) 2 /(gE^ lp ). They are thus small compared to the ones of the same structure in Eqs. ([3]) and ([7]) by the factor 
~ (Ej) 2 /[Ec i?™ p (l — a/%/2)] and can be neglected, unless one is in the vicinity of the isotropic Kondo model, 
a 2 «2. 



SM.B. EQUIVALENCE OF DIFFERENT FORMULATIONS OF ELASTIC SCATTERING 



In this section we address the relation between the current transmission amplitude, Eq. (|10|) in the main text, and 
the single photon elastic coefficients. The latter are encapsulated in the behavior of the time-ordered single-photon 
Green function Gt'\t{x'\x;u}) (with £,£' — L,R), where Gt'\i{x'\x;t) = —i{Tp{x',t)p(x,0)), 7 being the time-ordering 
operator [S4| . This Green function is related to the corresponding propagator 5u'G M (x'\x; u>) for a half-infinite lead 
detached from the impurity (Ej' : = 0, Cl,r = 0) by 

9 

Ge\i(x'\x;u) = 6 u >G h \x'\x;u>) -G U (x'\0;u}) T^Uco)G h \0\x; u), (S18) 

guj 

where the elastic T-matrix Tf}u{oS) has the structure given by Eq. © in the main text. The elastic T-matrix appears 

with a prefactor irv 2 /(gco) in Eq. (|S 18[) to compensate for the prefactors in the expansion of pe(x) in terms of the 
photon creation and annihilation operators [eigenmodes of the lead Hamiltonian, Eq. @ in the main text] , which 
reads (for a lead of length L with no-current boundary condition, d x (f>s(0) = 0, when decoupled from the impurity) 

Ps{x)= ^ COS (l X ) { a s,q - at,q) ( S19 ) 

q= 7rn/L. 

At zero temperature time-ordered and retarded Green functions arc the same for positive frequencies, uj > 0. There- 
fore, Eq. (|S18|) yields for the transmission coefficient ii(cj) for photons coming from the left, 

tL{u) = 2m ^L (PR^PlML (S20) 

™ 2 (PR(x%PR(0)t ((PL(0),p L (x)))™ 

Now, for the lead Hamiltonian, Eq. © in the main text, one has [S5| 



( Pe (x),pe(0))) 



hi = ^ iux/v 



irv 2 



(pR(x R ),p L (xL)t } =^Le<"l**-*>V«, (S21) 

where the superscript (0) denotes the propagator for a uniform waveguide, with no quantum impurities, as in the 
main text. Thus, Eq. (|S20I) can be rewritten as: 

. / \ ((pr{x'),Pl{x))) u {d x (j) R (x'),p L (x)))„ 
(Pr(x') iPl (x)))^ {dMx'),PL(x)))^ 
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where the last equality results from the equation of motion dtpe(x,t) — (vg / -K)d 2 4)i(x,t). We have thus proven the 
equivalence of Eq. (fTUf in the main text with Eq. (|S 18[) at zero temperature. The former can thus be thought of 
as a generalization of the latter to nonzero temperatures. Similar treatment applies to the other elastic scattering 
coefficients. 



SM.C. SHIBA RELATIONS FROM PHOTON SCATTERING 



The Shiba relation connects the low frequency behavior of the real and imaginary part of the Kondo local spin 
susceptibility Xzz(^) [sH IS7l. [S3l] . In this Section we will show how considerations based on photon-scattering can be 
used to rederive, as well as to generalize, this relation. Expanding the local spin susceptibility in powers of cj, 

He[x«(w)] =Xo + X20J 2 + ■ ■ ■ , (S23) 

Im[ Xzz (u)] =xiu + X3^ 3 + ■ ■ ■ , (S24) 

and substituting in Eq. (j!4|) in the main text, we obtain an expansion of the total inelastic scattering probability 
ji(u>) in powers of u <C Tr-, with coefficients depending on the %%■ O n the other hand, from Eq. (|18[) and Eq. (|15[) in 
the main text it follows that when u> is small, 7^(w) ~ oj 4 in the presence of a magnetic field, while r yi(u)) ~ w 6 for 
B z = 0. Comparing these results we find that the vanishing of the total inelastic scattering probability je(uj) to order 
uj 2 leads to the Shiba relation [cf. Eq. (|T2|) in the main text] [S7j : 

Xi=7ra 2 ( Xo ) 2 , (S25) 

whereas the vanishing of 7^(w) to order w 4 in the absence of a magnetic field leads to a new, higher order, Shiba-like 
relation: 

X3 = ira 2 [2xoX2 + (xi) 2 ]- (S26) 
This latter relation can be easily verified to hold at the exactly-solvable Toulouse point a = 1 [H, [13, IH| . 



SM.D. INELASTIC SPECTRUM FROM NONLINEAR RESPONSE FUNCTIONS 



The inelastic spectrum 7^u(u/|a;)du/ is defined in the main text as the average number of photons within a frequency 
interval do/ around uj 1 emitted through lead £' for each incoming photons at frequency uj in lead £. Thus, it is a sum 
over the partial cross sections for all the possible multiphoton scattering processes, integrated over all the photons 
except the one with frequency ui'. In this Section we will show that, similarly to Eqs. (|10l) - (|ll[) in the main text for 
elastic scattering, 7^ (u/|u;) can also be expressed in terms of response functions, and related to local spin correlators. 
Since the number of photons emitted at frequency ui 1 is proportional to the flux of incoming photons, or incoming 
energy flux (assuming scattering between two or more incoming photons is negligible), we need to consider second 
order response to the ac source of incoming photons. 

To spare us the need to carry around the indices £, £' in the following, we define 

7(wV) = E >K(^V), 7/'|<(<"V) = ^T7(wV), (S27) 
— ' or 
l,V 

The second relation results from the fact that only the "spin fields" 4> s (x) and p s (x) are coupled to the impurity [cf. 
the discussion following Eq. ^ in the main text]. 

The quantity of interest here is the time-averaged rate of change of the photon number n s ^ — a\ k ,a s ,k> , k' = co'/v, 
to second order in an applied ac voltage H ac = V{t)e 1,t p s {xi n ), with V(t) = 2Vbcos(u;£), r] —> + , divided by the 
incoming flux of photons of frequency uj [S8]. The photon creation and annihilation operators a\ q and a s ^ q are the 
Fourier modes of the bosonic fields (obeying the no-current boundary condition d x (f> s (0) = when decoupled from 
the impurity), 



q=7rn/ L, 
tj>0 



(x) = E i \[~^E C0S ( qx ^ ^ s ' q ~ a *>«) 



(S28) 



(S29) 



q—im/L 
q>0 
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where L is the lead length. The second order Kubo formula reads 



(n s ,k'(t)) 



(2) 



dt' 



diW ^A(. 4 .)ft(^.)(*-«'.*-*")n*')n*"). 



(S30) 



where the second order response function is: 

G ZlnM^yM Xi n) (* - f - *") = ~*(* - - *") ([[ ft -.k' (*)' P>( x ™ *')] . P.(*in, OD + {*' <* • (S3i) 

Here c, g denote "classical" and "quantum" fields in the Keldysh formalism [S9l. ISlClj ]. i.e., the sum and difference, 
respectively, of fields on the forward and backward contours. The time-average of the photon production rate can 
thus be written as: 



{n s ,k' 



v(2) 



dt 



(S32) 



where frequency arguments are in correspondence with time arguments in the previous equation. Since we should take 
the limit 77 — > + , the factor of 77 in this formula implies that we should be looking for contributions to the correlation 
function which are singular in that limit. Multiplying d(n Si fc/)( 2 ) / 'dt by the photon density of states L/(irv), and 
dividing by w|Vb| 2 / {nv 2 ), the rate of creation of photons with frequency uj moving towards the impurity by the source 
V(t), we have 



2r)vL 



- 1 \ - 1 \(w + w?, — oj + iri), 



(S33) 



For subsequent calculations it is better to look at a more general correlation function, where a\ k , and a Sy k' have 
different time arguments t\ and t2, respectively, and take the limit of t\ = 7)2 = t onl y at t he end. A suitable correlator 
is the following Green function, which appears naturally in the Keldysh formalism [SIO]: 



tl, t-i 



t\t 2 -t") = 

i9(h - t 2 )9(t 2 - - t") fe)} ,p s (x m ,t") 

+ ifl(t! - t')9{t' - t 2 )6{t 2 - t") ( [{ [aj ifc ,(*i),p a (a; in ,t / )] - ^^2)} , p s (x in) t")] ) 

+ i6(h -t')6{t' - t")8{t" - *2) ( { [ [a{ k , {h ), /5 S (x m , O] , Ps (x in , *")] . <k> (*a) } ) 

+ «i 2 ,/et"}. (S34) 

This formula reveals the general structure of the Keldysh Green functions for bosonic operators: all time orderings are 
allowed, provided the leftmost operator is classical, and each classical (quantum) operator appears in a commutator 
(anticommutator) with the operators to its left, i.e., the operators with larger time arguments. This structure will 
become important in the perturbative calculations in the next Section. Thus we can write 



y{u)'\w) = 



irjvL 



2tt 



(17, uj + irj 7 —uj + 777), 



(S35) 



where again frequency arguments are in correspondence with time arguments in the previous equation (the same 
convention will be followed for other four-point functions below) . 

Since the unitary transformations [such as li, defined before Eq. ([5]) in the main text] used to map between the 
different forms of the quantum impurity Hamiltonian [Eqs. ©, ©, and (J7J in the main text, as well as Eq. (|S49I) ] 
change the values of the charge densities pi(x) or the current densities oc d x 4>e(x) only locally, at x = 0, they do 
not affect the definitions of the scattering amplitudes. Thus, the inelastic spectrum 7(0/ \u>) [defined by Eq. (|S33j) or 
Eq. (|S35j) ]. as well as the elastic scattering amplitudes rz{uS) and ti(w) [defined by Eq. (fT0|) in the main text], can be 
calculated using any of these forms of the Hamiltonian which is more convenient. In the rest of this Section we will 
employ the spin-boson Hamiltonian, Eq. ([5]) in the main text. 

We will now show how the four-point Keldysh correlator appearing in Eq. (|S35[) can be written in terms of local 
four-point spin correlation functions. This can be done using the Keldysh path integral formalism, and integrating out 
the lead degrees of freedom. Alternatively, one may apply Keldysh perturbation theory to all orders in the spin-boson 
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FIG. SI. Schematic representation of the diagrams contributing to the four-point Keldysh Green function (|S34|) . Lines 
represent decoupled-lead propagators, and filled circles correspond to the spin-boson coupling H a = —7vvap a (0)S z [cf. Eq. (5}]. 
The dashed ellipses and square represent two and four point correlation functions of S z , respectively, calculated to all orders 
in H a . See the text for further details. 



coupling term, H a = —irvap s (Q)Sz [Sllj . Then, the four-point correlator defined by Eq. (|S34j) can be written as a 
sum of disconnected and connected diagrams, as depicted in Fig. IS 1 1 The former represent elastic scattering, and 
therefore vanish unless lj' = to. They can be shown to reproduce the square of the absolute values of the elastic 
scattering coefficients, Eqs. ([9|- (fTTj) in the main text. Since we are concerned here with inelastic scattering, we will 
rather concentrate only on the connected diagrams. These can be written as a product of four legs, representing two- 
point correlation functions of the lead operators calculated for a decoupled lead, multiplied by a four-point connected 
correlation function of S z , calculated with the full spin-boson Hamiltonian, Eq. ([5]) in the main text. We thus arrive 
at the following expression: 

Q ccqq + iii, —lj + iri) = ir 4 v 4 a 4 G hl '? c > - (— to — Z77)G- 1 '? C 1 - , n ,(oj — iv) 

a[ l? ,ia s . k r,p s (x irl );p 3 (x in ) y ' " /; Ps (x in );p s (0) V U p B (x in )-p s (0) " 

G 5t' C J;p a (0)^ + U + " ? ) G ? C J;P s (0)( -fi ~ W + i ^ G s7-k;S,;sS^ UJ + ir ), ~W + *V) 
+ G l\,;p 3 (0)( n + UJ + ^iZ^iO)^ - U + iv l) GC ™S z ;S^sS n ^ + ir >, -W + tV) 

+G™r. M0) iV + " + i V )G^. M0) (-n -lo + iv)&sSs. iS .;S M ( Q > w + ir,, -co + ir,)} . (S36) 

Since the four-point spin correlators are automatically connected for to' = to, we do not need to specify this explicitly. 
The disconnected lead correlation appearing in the above equation are given by 

G «£) ; «o,(")=S e '°""'"' (S37) 



G.t ° c _ , , (00) = coth —— 



(S40) 



where the last equation stems from the fluctuation-dissipation theorem. A similar relation holds for G^' cc .~ , Q s(co). 

The only terms in Eq. (|S36j) that are singular, and thus survive when the limit r\ — > + is taken in Eq. (|S35|) , are 
those which contain the product G^ cq (Cl + to + ir))G~' c ^.p ^(—Cl — to + ir)). By Eqs. (|S38[) — ()S39() . this product, 



together with the prefactor of 77 from Eq. (|S35[) . gives in that limit 

¥%r. A(0) (n + " + ir,)Gf^.. B{Q) {-n -u> + i v )^ ^-s(n + u + J). (S4i) 

One may then immediately perform the integral over Q. in Eq. (fS35)l . Plugging the result into Eq. (|S27| we are left 
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FIG. S2. = £\ e , =L R 7^/|^(w'|w), representing the inelastic spectrum summed over the incoming and outgoing leads, 

normalized by the total probability summed over the incoming lead, "/(to) = X^=l fi7 f ( aj )> a ^ the Toulouse point a = 1, for 
uj/Tk = 10.0, and B z — T = 0. The continuous line is the exact result for 7(a/|u;)/7(u;). The numerator here follows from 
Eqs. (|S42|I and (|S47|I evaluated numerically; the denominator is evaluated numerically with the help of Eqs. (|S46|I and Eq. (|14p 
in the main text. The dashed line corresponds to Eq. (|17|) in the main text, valid for uj' ,ui — uj' 3> Tk- Note the monotonic 
behavior for uj' > Tk, as compared with Fig. IS3I See the text for further details. 



with: 



Ji>\i(u'\uj) = —ala%u)u'x 



G S^;S 2 ;S i ( a; ' + W '- W ' W )- COth ( 



QCqqq 



(a/ + cj, — w, uj) — G q s cqq s . s . s (uj' + tjj, —uj, uj) 



(S42) 



Thus, inelastic scattering involves higher order local response function correlators than the elastic amplitudes: 
G s cqq s . s . s , the second order response of S z -S z correlations to the application of a local magnetic field, as well 
as G s qqq s , s . s , the third order local spin susceptibility. It can thus yield more information about the quantum 
impurity dynamics than elastic scattering can. 

The r.h.s. of Eq. (|S42|) can be evaluated exactly at the Toulouse point, a = 1, where the Kondo problem [ Eq. ([7]) 
in the main text] is equivalent to a noninteracting resonant level coupled to a spinless fermionic bath [S5l |S7| , [S3j , 



-Hrlm = v kc\c k + e (d)d - ±) + t Q d) ^ c fc + H.c, 

k 



k 



(S43) 



where eft (c£) creates an electron in the resonant level (mode k of the bath), with S z — > Sd—1/2 (and thus Eq = —B z ), 
as well as to = Ixy/ (2v / 2-7ra). The level width is V = tl/(2v) — I% y / (167rav). Since this model is quadratic, correlation 
functions of S z are easily calculated, using Wick's theorem and the results 

1 



G d:S (W) 



UJ ■ 



e ± iT 



tanh ( — ) 
\2TJ 



(S44) 
(S45) 



for the retarded, advanced, and Keldysh fermionic level Green functions, respectively [S£ 
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The dynamic spin susceptibility is then given by 

oo 

dn r 



2tt [{si + u- e ) 2 + r 2 ] [(n - e ) 2 + r 2 ] 



(Q + uj — eo — iT) tanh 



+ w 
2T 



(<? - _-„ + /Dial).]] ( — 



(S46) 

At zero temperature and magnetic field we get for the static susceptibility, Xzz(O) = l/(7iT). Thus, Tk = 7iT according 
to our definition [cf. the discussion before Eq. of the main text]. In addition, the correlation functions appearing 
in Eq. (|S42|) are given by 



G 



abqq 

S z ]S z'-S z';S ; 



(uj + uj' , —u), uj) = —i 



2^ 



-Tr 



with a,b = c, q, and where 



G(u>) = 



f a G(uj' + n)T b G(n)G{uj + O)G(fi) 

+ T a G{uj' + Cl)G(u + uj' + VL)T b G(uj + n)G(n) 
+ T a G(uj' + Q)G(oj + uj' + fl)G(uj' + Q)f h G(n) 
+ {uj < Uj} , 







(S47) 



(S48) 



and f c is the Pauli matrix t x , whereas f 9 is the unit matrix. An example of the resulting inelastic spectrum is plotted 
in Fig. E2 



SM.E. PERTURB ATIVE CALCULATION OF THE INELASTIC SPECTRUM 



When u), uj' 1 and \uj — uj'\ or T are large with respect to the Kondo temperature, one may evaluate the inelastic 
spectrum j£i\g(uj'\ui) perturbatively in I xy oc EjEf [cf. Eq. (|4]) in the main text] for any value of a, and obtain 
Eq. (|16[) in the main text. In this section we will present the details of this calculation. 

In this regime it is useful to apply the transformation H — » V^HV with V = e la ^ s <°^ s * to the spin-boson Hamiltonian, 
Eq. ([5]) in the main text (as argued in the previous Section, such a transformation does not affect the scattering 
amplitudes), so as to transfer the impurity- leads coupling into the perturbative I xy term, 



(S49) 



H * = E ^ / { [ d M x )} + kPA^)] 2 } dx - B Z S Z - ( e -**-<°>S + + e **(°)£L 



Expanding the Keldysh Green functions appearing in Eq. (|S42I) in I xy , the zeroth and first order terms vanish. 
The second order terms breaks down into a products two-point boson correlator and a six-point spin correlator, to be 
evaluated for the Hamiltonian (|S49|) with I xy = 0: 



QCcqq 



(t 1 -t 2 ,t 1 -t\t 2 -t")=i(J^ X 



a,b—q,c 



(S50) 



]T / d*i / d aa (^(t 1 )^(t a )5»(O^(t , 0^(*i)^(* 2 )>([e- i ^'( o '' 1 )] B [ e fa *-( o '' ! ')] S ) ) 



oo — oo 



Qcqqq 



(S51) 



£ J d Sl I d S2 (^(t 1 )5K^)^KiO^(i'')^( S i)^( S2 ))([e-^( ^)] a [e^^)] 5 ), 

a,b—q,c 



-oo — oc 



whereas G q g qq s . s . s (t\ — t 2 ,ti — t',t 2 — t") is obtained from G c ™ s . s . s (t 2 — ti,t 2 — t',ti — t") by interchanging t\ 
and t 2 . Here a — q,c for a = c, <j, respectively, and similarly for b. Therefore, the term with a = b = c contains a q-q 
lead correlator, and thus vanishes. 
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For oj, oj', \oj — oj'\ 3> B z , one may neglect the effects of the magnetic field. Then, th e spi n operators appearing in 
Eqs. (|S50[) - (|S51[) are time independent. Following the rules of the Keldysh formalism |S10| . the corresponding spin 
correlators can be written as combinations of commutators and anticommutators of the spin operators, depending 
on the ordering of the time arguments. Most of these turn out to be zero. The spin correlator on the r.h.s. of the 
Eq. (|S50|) does not vanish only if a = b = q, in which case it gives 

(Sl(ti)Sl{t2)SW)Sltf')Sl{ai)&_(82)) = 9{t 1 -s 1 )9{s 1 -t')9(s 1 -t")9(t'-s 2 )9(t"-s 2 )9(s 2 -t 2 )+{s 1 o s 2 }, (S52) 

whereas the spin correlator on the r.h.s. of the Eq. (|S51[) docs not vanish only if a = c, b = q, when 

(S c z (h)SUt2)S q z (t')S q z (t")S c + ( Sl )Si( S2 )) =9(h - s 2 )9(s 2 - t 2 )9{s 2 - t')9(s 2 - t")9(t 2 - Sl )9(t' - Sl )9{t" - 8l )+ 

9( Sl - t 2 )9( Sl - t')9{ Sl - t")9{t 2 - s 2 )9{t' - s 2 )9(t" - s 2 )9(s 2 - t x ), 

(S53) 

or if a = q, b = c, in which case one should simply interchange si and s 2 in the last equation. 

Plugging Eqs. (|S52[) - (|S53|) back into Eqs. (|S50p - (|S51|) . one can perform the integrals over S\ and s 2 , and then 
calculate the Fourier-transform of the results. Using in addition the fluctuation-dissipation theorem to express all the 

different Keldysh lead correlators in terms of the retarded one, = (le la ^ 3 ^; e~ lQ * s (°) » , we find 



GZ% M (" + =^S { 2c ° th (|Q Im - C ° th i^f) Im 

(S54) 

G %%s mi s n -,s m ( u + w',-w,w) = - G q s c ™ Sz . s ^. s ^ + w'.-w.w) = 

JjS (fe) 2 { 2Im Im + ^ + Im - ^ I' (S55) 



Substituting these expressions into Eq. (IS42[) we arrive at Eq. (IT5|) in the main text: 
\47ra / 



7/'|i(wV) - (^) [im - «')] - + «b(w')] [1 + n B{u - oj')] - BB^n^u - a/)} 

+0(u/ - w)|n B (a;') [1 + n B (w' - u>)] - [1 + n B {oj')] n B (oj' - 
+Im [x+_(w + w')] |n s (w + a;') [1 + n B (cj')] - [1 + n B {w + oj')] n B (oj')\ , 



where [SI 



hi , x 1 f™ 2 \ 1 /2ttT\ q ~ x fa 2 fi 



with 5(x,y) the Euler beta function [S12] . Thus, Im[xf_(tt)] oc [max(fi, T))] Q - 1 for small and T. As mentioned 
in the main text, the first line of Eq. (fl~6| describes a process where an incoming photon at frequency u is absorbed by 
the quantum impurity, and a photon at frequency oj' < oj, plus additional photons whose energies sum up to oj — oj' 
are emitted (this is the only process allowed at zero temperature) and the reverse process. Similarly, the second line 
describes a process where a photon at frequency oj' > oj is absorbed, and a photon at frequency oj, as well as photons 
whose frequencies sum up to oj' — uj are emitted and vice versa. Finally, the third line describes a process where 
photons whose frequencies sum up to oj' + lj are absorbed, and photons at frequencies oj and oj' are emitted and vice 
versa. 

The structure of Eq. (|16l) in the main text thus suggests that it can be obtained from a kinetic equation. Indeed, 



one can write down the Boltzmann equation for the average mode occupations n q = (a| g a Sjg ), accounting for all the 
possible multiphoton scattering processes to second order in I xy . The corresponding probabilities can be obtained by 
Fermi's golden rule from the Hamiltonian (IS49[) . after expanding the exponents in the last term of the Hamiltonian 
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to all orders in the bosonic fields |S13| |: 




where ljq = vQ. In equilibrium [in the absence of the external driving V(t)], the mode occupations are given by the 
Bose-Einstein distribution, uq = n B (ujQ). In order to find the rate of change of occupation of mode q — k' by the 
ac excitation V{t) to second order in I xy , one should substitute on the right hand side of Eq. (|S57|) the equilibrium 
(Bose-Einstein) occupations for all the modes, except for the mode with wavevector k — uj/v, whose occupation is 
modified by n% by the ac source V(t); thus, uq — ti b {ujq) + {it / L)n\5{Q — k). Multiplying the resulting rate by 
the photon density of states L/(irv), and dividing by the incoming flux of photons of frequency u [i.e., n%v/ (2L)], we 
recover our previous result, Eq. (1161) in the main text, if we employ the following relations: 



[1 + n fl (£))] Im^e <Q ^°);e~^ (0) )) = * 



a N+N' 



N,N'=1 



n\n'\ 



dgi 



dgjy f dq[ 
QN J q[ 



<i'n> 



n B {w qi ) ■ ■ ■ n B {u qN ) [l + n B (u q ^)\ ■•• [l + n B (u q ' N ,) 6{tt + uj qi H h 



Uq N -L) -U 



1l 



-oyj, (S58) 



n B (fi)Im((e' 



hi 



/, ! ^(O). e -ia0(O)^\ 



oo 



a 



N+N' 



^ N\N'\ 

N,N' = 1 



dgi 
9i 



dq N f dq[ 



QN 



IN' 



n B (w qi )---n B (u qN )[l + n B (u q , i )]--- l + n B {u q , ) 5{w qi H h 



- u„„ - u) - u)„ 



- ay - 0). (S59) 



SM.F. INELASTIC SCATTERING IN THE SMALL a LIMIT 



In this Section we will analyze inelastic photon scattering in the limit of small a at zero temperature. In that 
regime it is useful to use the spin-boson version of the Hamiltonian, Eq. in the main text. We will start from 
the case B z = 0. Then we have a two-level system (S x = ±1/2), where the two levels are separated by Tk — Ej, 
and weakly coupled to the bath of photons. Since every photon emission or absorption flips the impurity spin, the 
inelastic process which is lowest-order in a and leaves the two-level system in its ground state involves four photons. 
The amplitude for a photon at frequency uj incoming in lead £ to scatter into photons of frequencies a/, u\, and 
LU2 = io — lo' — uj\ outgoing into leads £' ', i\ and £2, respectively, is, to the lowest order in a, a sum over the partial 
amplitudes of the 4! = 24 different orderings of the absorption of the single incoming photon and the emission of the 
three outgoing ones. Squaring this total amplitude and multiplying by the appropriate density of states factors we 
find the cross section 



UJUJ UJ1UJ2 



(w - E^ r ){uj + E^ R )(w' - E L j R )(lu> + E^ R )(w! - £^)(u;i + E L j R ){uj 2 - E LR )(uj2 + E$ R ) 



(S60) 



where Ej R 



Ej R [l — (a 2 /2) \n(bj/Ej R )\ +iTj R accounts for the shift and finite lifetime (broadening) of the excited 



impurity state, with Tj R = ira 2 Ej R /4 (The shift in the real part of Ej R corresponds to the change in Tk, Eq. 
in the main text, calculated to order a 2 ). Thus, while the cross section is only of order a 8 for small a (for fixed 
ai/aji), it displays peaks of height cx a 4 and width cx a 2 whenever one of the frequencies is close to Ej R . It should 
be noted that having more than one of the outgoing frequencies u/, o>i, and LU2 close to Ej R does not lead to even 
higher peaks, since the numerator in Eq. (|S60|) vanishes in that case. 
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Eq. (S61) 

Eq. (S62) 

Eq.(18) 




FIG. S3. y(u)'\uj) = £\ e , =L R j£ii£(oj'\ui), representing the inelastic spectrum summed over the incoming and outgoing leads, 
normalized by the total probability summed over the incoming lead, j(u) = Ef =i .i{7^(^), for « 2 = 0.1, uj/E^ r = 30.0, and 
B z = T = 0. In this regime 7(0;) ~ a e '(Ej R /uj) 2 \n(ui/Ej R ). The inset is a zoom-in into the region u) - J < E LR . The 
continuous line is the exact result (to leading order in a), obtained from numerical evaluation of Eq. (|S61[) together with 
Eq. ()S60[1 . The dashed line corresponds to Eq. (|S62|) . valid for Ej R <C U)' < uj — Ej R , and the dotted line corresponds to 
Eq. (|19[) in the main text, valid for lu — uj' <C Ej . Note the nonmonotonic behavior for uj' > Tk, as compared with Fig. IS2I 
See the text for further details. 



Integrating over U3% and summing over tx t 2 we find the four-photon process contribution to the inelastic spectrum 

7<^(w'|w) = ^2 / T^'AAI^'j^i.^ - u ' ~ wi|w)da;i. (S61) 

t u ta=L,R 

Let us discuss the main features in the dependence of j^ e (uj'\uj) on uj, uj' and a. When all frequencies are small 

with respect to E^ R , no resonance contributes, leading to the second term in Eq. (fT5)l in the main text, with 
a u (a) = 37r 2 a 4 /4. In that case, therefore, the spectrum 7^ (uj'\uj), as well as the total inelastic scattering probability 
7^(oj) [obtained through the sum rule, Eq. (|15p in the main text], are very small, of order a s . 

For uj ^> Ej R (more precisely, uj > 2Ej R ), the behavior is richer, as depicted in Fig. IS3I For uj' > uj — Ej R none 
of the frequencies is close to a pole, and the spectrum is still cx a 8 , corresponding to the second term in Eq. (|19| in 
the main text, with a' u {a) = ir 2 a /12. For uj' < uj - Ej R the integration over uj% includes the regions 10^2 « Ej, so 
the spectrum is oc a 6 in most of this range, except for a peak of height cx a 4 and width oc a 2 when uj' ~ Ej R . Away 
from that peak, in the regime Ej <C uj' < uj — Ej R the calculation can be carried out explicitly to the lowest order 
in a, leading to, 

(4)/ /I \ 2 2 2 2 ( E J R ) — Uj' — Ej R ) 

UJUJ (uj - UJ') 

It should be noted that Eq. (|S62j) agrees with Eq. (fl7|) in the main text in their common domain of applicability, i.e., 
lowest order in a and the range Ej R <C uj' <C uj — Ej R . By Eq. (fT5|) in the main text, this latter range gives the 
dominant contribution to the total inelastic scattering probability je(uj) for u ~> Ej R . Eq. (|S62[) results in 7^(w) ~ 
(Ej R /ui) 2 a 6 \n(uj/E 1 j R ) for a 2 ln(ui/Ej R ) < 1, whereas Eq. (J7]) in the main text shows that ^(w) ~ a 4 (Ej R /uj) 2 
for a 2 ln{uj/Ej R ) > 1. 

In the regime Ej R < uj < 2Ej R a similar analysis leads to a total inelastic probability cx a 6 . Finally, when uj itself 
is resonant, uj sa Ej R , both 'ye'\t(ui'\u>) and 7f (w) are cx a 4 . The peaks when one of the frequencies uj', uji, uji, is also 
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close to Ej R are suppressed here by the frequency factors in the first line of Eq. (|S60I) . since the other two frequencies 
must be close to zero in this case. 

Turning on a finite magnetic field B z , three-photon processes become possible. To lowest order in B z /Ej R their 
contribution to the inelastic spectrum is 



7^(cj'|cj) = '—a\a.\,a. 



2„2 2 D 2 



Bi 



E 



LR 



wu'(w — CJ ) X 



cj 2 + cj' 2 - cjcj' - 3 (Ej R 



(cj - E$ R ){u + E^ r )(oj' - E L j R )(lu' + E$ R )(w - cj' - E L j R ){lu - cj' + Ej R ) 



(S63) 



At small frequencies we now recover the first term in Eq. (|18|) in the main text, with as (a) = 97r 2 a 2 /2, i.e., 7£<^(cj'|cj) oc 
a 6 . For cj > Ej R the spectrum has two peaks, at cj' sa Ej R and cj — cj' « Sj*, both of height and width oc a 2 , 
leading to total inelastic probability oc a 4 , whereas for cj — cj' <C we recover the first term in Eq. (fT9|) in the 
main text for cj 3> Ej R , with o' B (a) = 7r 2 a 2 /2. Finally, for cj -Ej^ we have 7^^(cj'|cj) oc a 2 , with narrow peaks at 
cj' ~ Tj R and cj — cj' ~ resulting in 7^(cj) oc a 2 ln(l/a 2 ). On the other hand, all these values are suppressed by 
a factor - (B z /Ej R ) 2 . 
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